Mechanical surface tension governs membrane thermal fluctuations 
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Motivated by the still ongoing debate about the various possible meanings of the term surface 
tension of bilayer membranes, we present here a detailed discussion that explains the differences be- 
tween the "intrinsic" , "renormalized" , and "mechanical" tensions. We use analytical considerations 
and computer simulations to show that the membrane spectrum of thermal fluctuations is governed 
by the mechanical and not the intrinsic tension. Our study highlights the fact that the commonly 
used quadratic approximation of Helfrich effective Hamiltonian is not rotationally invariant. We 
demonstrate that this non-physical feature leads to a calculated mechanical tension that differs dra- 
matically from the correct mechanical tension. Specifically, our results suggest that the mechanical 
and intrinsic tensions vanish simultaneously, which contradicts recent theoretical predictions derived 
for the approximated Hamiltonian. 
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I. INTRODUCTION 

Bilayer membranes are quasi-two-dimensional (2D) 
fluid sheets formed by spontaneous self-assembly of lipid 
molecules in water Their elasticity is traditionally 
studied in the framework of the Helfrich effective surface 
Hamiltonian for 2D manifolds with local principle curva- 
tures ci and C2 [2J 



Ho = / dS 



1 , n2 

co + (ci + c 2 - 2c ) +k cic 2 



(1) 



where the integration is carried over the whole surface of 
the membrane. The Helfrich Hamiltonian involves four 
parameters: the spontaneous curvature cq, the surface 
tension cro, the bending modulus koj an d the saddle-splay 
modulus Rq. For symmetric bilayer membranes, Co = 0. 
If, in addition, the discussion is limited to deformations 
that preserve the topology of the membrane, then (by 
virtue of the Gauss-Bonnet theorem) the total energy as- 
sociated with the last term is a constant, and one arrives 
to the more simple form 



Ho= / dS 

J A 



1 2 

0-0 + (ci + c 2 y 



= <j A+Ik J 2 , (2) 



where A is the total area of the membrane and J, defined 
by J 2 = J dS (ci + C2) 2 , is the integrated total curvature. 

The surface tension appearing in Eq.Q is known as 
the "intrinsic tension". It represents the elastic energy 
required to increase the surface of the membrane by a 
unit area, and can be identified with the derivative of 
the energy with respect to the area A, <r — (dH/dA), at 
constant J. As this quantity is not directly measurable, 
its physical meaning is still a matter of a fierce debate. In 
molecular simulations, one can attach the membrane to a 
"frame" and measure the "mechanical (frame) tension" , 
r, which is the lateral force per unit length exerted on the 
boundaries of the membrane . Experimentally, the 
mechanical tension is routinely measured by micropipette 
aspiration of vesicles @, Formally, the mechanical 



tension is obtained by taking the full derivative of the 
free energy F with respect to the frame (projected) area: 
r = dF/dA p @. From a comparison of the above def- 
initions of do and r, it becomes clear that the intrinsic 
and mechanical tensions are different quantities. Con- 
trary the former, the latter is a thermodynamic quantity 
that also depends on the entropy of the membrane. Bi- 
layer membranes usually exhibit relatively large thermal 
undulations at room temperature [lol | and, indeed, their 
mechanical tension also includes an entropic contribution 
due to the suppression of the amplitude of the undula- 
tions upon increasing the projected area. 

The surface tension can be also measured indirectly by 
recording and analyzing the statistics of the membrane 
height fluctuations Ullll. The analysis is based on the 



so called Monge parametrization, where the surface of the 
fluctuating membrane is represented by a height func- 
tion, h(x,y), above the frame (x,y) plane. The Helfrich 
Hamiltonian does not have a simple form when expressed 
in terms of h(x, y). However, for a nearly flat membrane, 
i.e., when the derivatives of h with respect to x and y 
are small (|c\;/i|, \d y h\ <C 1), one obtains the quadratic 
approximation 
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Note that unlike Eq.©, the integral in Eq.Q runs over 
the frame area rather than over the area of the man- 
ifold. The quadratic approximation can be diagonal- 
ized by introducing the Fourier transformation: h q = 
(1/Ap) J dxdy h(x,y) exp(—iq ■ r). In Fourier space the 
Hamiltonian reads 

At, 
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n q ) \h q \ 



(4) 



and by invoking the equipartition theorem, we find that 
the mean square amplitude of mode q ( "spectral inten- 
sity"): 

k B T 



(\h q \ 2 ) 



A p (a q 2 + n q 4 ) ' 



(5) 
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From the last result, it seems as if <7o can be extracted 
from the fluctuation spectrum of the membrane. Re- 
sults of both fully- atomistic and coarse-grained simula- 
tion (l3l - [T7t show that spectral intensity can indeed be 
fitted to the form 



(Kl 



A p (rq*+0(q*)) 



(6) 



However, the derivation of Eq.© is based on the ap- 
proximated Hamiltonian Wf and, therefore, it is not a- 
priory clear why the so called "^-coefficient" appearing 
in Eq. (|6]) , r = Co- in fact, some theoretical studies have 
argued that the q 2 coefficient is actually equal to the 
mechanical tension r [I8I - I20I ] . This conclusion has been 
rejected more recently in favor of the more common in- 
terpretation that r = <t [HI HH • Membrane simulations 
carried at fixed mechanical tension (usually performed 
for r = 0)tend to agree with the result that r = r 
d, HH, [H, [H, [23| , but simulations that show the oppo- 
site r 7^ r also exist [U [24|. In this paper we settle 
this problem and prove, analytically and computation- 
ally, that the correct result is indeed r = r. 



II. WHAT DOES THE INTRINSIC TENSION 
REPRESENT? 

Much of the confusion associated with the physical 
meaning of surface tension in membranes is related to 
fact that the concept of surface tension has been origi- 
nally defined for an interface between bulk phases (e.g., 
between water and oil) [25[. In its original context, the 
surface tension represents the access free energy AF per 
unit area A of the interface between the bulk phases 
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AF 



(7) 



When the concept of surface tension is introduced into 
the theory of bilayer membranes, its meaning is distorted 
due to the following two major differences between mem- 
branes and interfaces of bulk phases 

1. In the case of bulk phases, it is often assumed that 
the interface between them is flat. There is usu- 
ally very little interest in the thermal roughness of 
the interface, unless it is very soft. In contrast, bi- 
layer membranes are treated as highly fluctuating 
surfaces whose elastic response is very much influ- 
enced by the entropy associated with the thermal 
fluctuations. 

2. The changes in the areas of both systems arise 
from very different origins. In the case of, say, 
a water-oil interface, the changes in the interfa- 
cial area are not produced by elastic deformations 
(dilation/compression) that modify the molecular 
densities of the bulk phases. Instead, they result 
from transfer of molecules between the bulk phases 



and the interface occurring, for instance, when the 
shape of the container is changed. The surface ten- 
sion 7 is essentially a chemical potential which is 
directly related to the exchange parameter between 
the coexisting phases [26| . The case of bilayer mem- 
branes is completely different. Here, there is an ex- 
change of water between the bulk fluid and interfa- 
cial region, but almost no exchange of lipid material 
because the concentration of free lipids in the em- 
bedding solution is extremely low (10 -6 — 1CP 10 M 
[l[). In other words, there is no reservoir of lipids 
outside of the bilayer and, therefore, a changes in 
the bilayer area result in a change in the area den- 
sity of the lipids. The membrane surface tension 
measures the response to this elastic deformation, 
including the indirect contribution due to the ex- 
change of water between the bilayer and solution 
resulting from the deformation (which has influence 
on the effective elastic moduli). 

Most earlier theoretical investigations of membranes 
involved the assumption that the area per lipid a is con- 
stant. This assumption relies on the observation that the 
energy cost involved in density fluctuations is much larger 
than the energy scale associated with curvature fluctua- 
tions. Further assuming that the lipids are insoluble in 
water (and, therefore, they all reside on the membrane 
where their number N is constant) implies that the to- 
tal area of the membrane A = Na is constant. If that 
is the case, then why does one need to include this con- 
stant in the Helfrich Hamiltonian [first term on the right 
hand side of Eq.©], and what does the coefficient, uo, 
represent? The answer is simple. It is technically very 
hard to calculate analytically the partition function for a 
fluctuating manifold with a fixed area A. The first term 
in Eq.© is a Lagrange multiplier that fixes the mean 
area (A) of the membrane, and the value of 00 is set by 
the requirement that A = (A) = dF/d<r n . As usual, it 
is assumed that in the thermodynamic limit the relative 
fluctuations in the total area become negligible, and there 
is no distinction between (A) and A. 

The renewed interest in the meaning of surface tension 
during the past decade is very much linked with the rapid 
development in computer modeling and simulations of 
bilayer membranes. In molecular simulations the num- 
ber of lipids N is usually fixed, but the total area in 
not. Membranes are no longer treated as incompressible 
thin films, but rather as stretchable/compressible sur- 
faces whose elastic response result from their intermolec- 
ular forces. With this point of view in mind, it is clear 
that the intrinsic tension cto in Helfrich Hamiltonian rep- 
resents an elastic coefficient which, potentially, may be 
related to the mechanical tension r. There is, of course, 
no reason to expect that the elastic energy E of the mem- 
brane is linear in A (27j . A quadratic elastic function 
E = 1/2Ka(A—Aq) 2 seems like a more appropriate form, 
where Ka is the stretching/compression modulus and Aq 
is the relaxed area of the membrane [28[ (also known as 
Schulman's area [Ij|). Expressing the total area as the 
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FIG. 1: A fluctuating membrane with a projected area A p 
at four different states: (a) the reference undeformed state, 
(b) under shear deformation, (c) subjected to an arbitrary 
deformation, and (d) as in (c) but with periodic boundary 
conditions. Each state is characterized by the mean height 
profile function h(x, y) which is indicated by a thick solid 
line. The dashed thin line in (a) serves as a reminder that the 
membrane is fluctuating around its mean profile. 



sum of the projected and undulations areas A = A P +SA, 
and assuming that 5 A <C | A p — Aq \ , yields the linear ap- 
proximation E ~ l/2K A (A p - A ) 2 + K A (A P - A )6A 
from which one identifies that <To = Ka{A p — Aq). The 
assumption that 5 A <C \A p — Aq\ becomes increasingly 
accurate at high tensions; but, nevertheless, the linear 
relationship between E and A is used in Eq.© over the 
entire range of intrinsic tensions. 



shear is a special case of a more general family of elastic 
deformations resulting from the action of normal stress 
fields. A representative of this family of deformations is 
shown schematically in fig. Q](c). Mathematically, such 
deformations can be characterized by the mean height 
profile of the membrane: (h(x,y)) = h(x,y). This func- 
tion, which vanishes in the undeformed reference state 
h(x, y) = 0, serves as the strain field that describes arbi- 
trary deformations of the fluctuating membrane at fixed 
projected area A p . Specifically for the shear deforma- 
tion shown in fig.[Hb), h(x,y) — ex where e is the shear 
strain. The normal stress field associated with the strain 
function_/i(a;, y) can be derived from the elastic free en- 
ergy F(h). Based on the same physical arguments used 
to justify the introduction of Hclfrich effective Hamilto- 
nian Eq.©, we hereby assume that the elastic free energy 
F (h) can be written in a similar form 
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F{h(x,y))=aA (h) + -k(J) 2 (h) , 



(8) 



where A(h) and J(h) are the total area and integrated 
total curvature of the mean height profile of the mem- 
brane, h(x,y). Our conjecture of Eq.® is based on 
the fact (which will be proved in the following section 
IIII Ap that this form of the free energy yields the experi- 
mentally and computationally well established Eq. © for 
the spectral intensity. The coefficient a and k appear- 
ing in Eq.Q are the renormalized surface tension and 
bending rigidity. These quantities are thermodynamic 
properties of the membranes which, in general, include 
entropic contributions and, therefore, are different than 
the corresponding microscopic coefficients appearing in 
the Hclfrich effective Hamiltonian. 



A. Linear response theory 



III. HELFRICH ELASTIC FREE ENERGY 

The proof that r = r consists of three steps and in- 
volves the introduction of two additional surface tensions: 
(i) (J, the renormalized surface tension, and (ii) //, the 
out-of-plane shear modulus (which, as it will turn out 
below, is also a surface tension). Let us first understand 
what these two quantities represent. The out-of-planc 
shear modulus [i is the force per unit length required to 
introduce the deformation depicted in fig. [Ub) from the 
initial reference state depicted in fig. [2(a) . This defor- 
mation is achieved by applying opposite normal forces 
on boundaries of the membrane. (Notice that the out- 
of-plane shear deformation changes the area of the mem- 
brane, but preserves the volume of the three-dimensional 
Euclidean metric. We will revisit this point later in sec- 
tion [nTC]) The shear modulus is not related to the de- 
formation of a specific configuration but rather to the de- 
formation of a fluctuating membrane. The out of plane 



The first step in the proof that r — t is to show 
that r = a, which follows from static linear response 
theory. Detailed proof of this point can be found in 
ref. [l9j ; here we provide a short version of the derivation. 
Let us consider an elastic deformation with a function 
h(x, y) that satisfies periodic boundary conditions (see, 
e.g., fig. [T|d)). Introducing h q , the Fourier transform of 
the function h(x,y), and expanding A and J in powers 
of h q , we find 

F ({h q }) =aA p + ^J2 + K ? 4 ] M-« + ° (N 4 ) 
<? 

(9) 

where = h* since the function h(x,y) is real [30]. 

Let us now consider the perturbed Hamiltonian % = 
Hq — "^Zq-iqhq, where Ho is the Hamiltonian of the mem- 
brane. Notice that we do not need to know the specific 
form of Ho and, in particular, we do not assume here 
that Ho is necessarily the Helfrich Hamiltonian Eq.([2]) or 
its quadratic approximation Eq.@. Instead, we assume 



that the free energy is given by Eq.© [or by its Fourier 
space counterpart Eq.©]; and as we shall now show, this 
is the key assumption that ultimately leads to Eq.(j6|). 
The derivation is as follows: Introducing the Gibbs free 
energy G({j q }) = -fc s Tln[Tr exp(-H/k B T)], and tak- 
ing the derivative of G with respect to j q yields 



h a 



dG({ Jq }) 

djq 



(10) 



The conjugate thermodynamic variables h q and j q can 
be also related to each other via 



dF({h q }) 



dh„ 



(11) 



which, by using Eq.© for the Hclmholtz free energy F, 
reads 



Jq 



A p (aq 2 



nq 4 ) h_ q 



0(\h q 



(12) 



Deriving Eq. (fT0|) with respect to j- q yields the relation- 
ship between the elastic linear response and the equilib- 
rium fluctuations in h„ 



(IM 2 ) 



i rrdhq 

KrT- 



dj-q 



j q =o 



dj-c, 
dh„ 



(13) 

where (• • • )o denote a thermal average using the unper- 
turbed Hamiltonian % . Using Eq. ([T2l) in Eq. (fT3)) . one 
arrives to the result that 



k B T 



A p (aq 2 + K q 4 )' 



(14) 



This completes the proof that r = a which, as noted 
above, is independent of the explicit form of "Ho- 



B. The renormalized tension is a shear modulus 

The second step is to prove that a — /i. This can 
be done easily by considering the shear deformation 
h(x,y) = ex shown in fig. Hlb). The out of plane shear 
modulus is defined by the expansion of the free energy 
density in powers of e 



F 



fo 



M 2 

2 £ 



0(e 4 ). 



(15) 



(The expansion includes only even powers of e because of 
the symmetry of the problem with respect to reflections 
in the direction normal to the projected area.). However, 
for the shear deformation, the mean profile is flat, i.e. 
J = 0, and therefore Eq.© reads 



F(h(x,y)) = uA. 



(16) 



The area A and e are related by A = vT + e 2 A p — A p [l- 
e 2 /2 + C(e 4 )], which yields 



F 



(17) 



unstrained 



sheared 



rotated 



FIG. 2: (a) The unstrained reference state of a thin film, 
(b) The film subjected to a simple shear deformation, (c) 
The film subjected to a pure shear deformation. The two 
deformed states have the same elastic free energy since they 
can be transformed into each other by rotation. The volume 
of the film (indicated by the gray shaded area) is the same at 
all three states. 



Comparison of Eqs. (p~5)) and (fl7|) leads to the result a = /j 
which, again, is independent of the explicit form of the 
membrane Hamiltonian "Ho- 



C. Rotational invariance 

Having demonstrated that r = a = fi, the last step 
in the proof that r = r is to show that r = a = [i. 
This follows immediately from the invariance of Hclfrich 
free energy Eq.([5]) with respect to rigid transformations. 
Specifically, the sheared membrane which is replottcd in 
fig-E^b), can be rotated (see fig. H[c)) so that the mean 
profile of the deformed membrane lies in the (x, y) plane 
defined by the undeformed state (fig. (Ha)). In this orien- 
tation A — A p , and from Eq. ([T6f one gets that the free 
energy cost of imposing the deformation in fig. [2jc) is 
F 

tension 



aA p . From this result one finds that the mechanical 



dF 
dA r , 



(18) 



While Eq. fTS)) is correct, it misses an important issue 
which cannot be captured within the framework of Hel- 
frich model that treats the membrane as a 2D manifold 
with no 3D volume. In figs.[2][a)-(c) we have intentionally 
drawn the membranes as thin films, which highlights a 
very important point. The transition from the reference 
state (a) to the deformed state (c) is achieved via a se- 
quence of two volume-preserving transformations - shear 
followed by rotation. Therefore, the reference state and 
the deformed one have different areas but the same vol- 
ume, which means that in Eq. (|18p one must take the 
derivative at a constant volume: 



dF 



(19) 



Experimentally and computationally, there is no prac- 
tical way to fix or even determine the volume of a bi- 
layer at the molecular resolution. This, however, does 
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not render the above argument irrelevant to lipid bilayer 
membranes. The membrane and the embedding fluid 
medium are placed in a "container" whose volume can be 
easily controlled. Eq. pl?)) states that the frame tension 
of the bilayer can be measured by changing the cross- 
sectional area of the container while keeping its volume 
fixed. This tension is associated with the entire "interfa- 
cial region" of the system that includes both the bilayer of 
lipids as well as the hydration layers of structured water 
around the bilayer. The fluid bulk water has no elastic 
response to volume preserving deformations and, there- 
fore, it make no contribution to r. 

The fact that r is associated with a volume preserving 
deformation is important. The reader may have noticed 
that in Eq. ([T8| we wrote r = a, while in Eq. (fT9l) we pre- 
ferred the equality r = fj,. Obviously, the latter is also 
correct since, as discussed above in section UlIBl er = \i. 
In writing r = (J, in Eq. (|19l) we wish to emphasize the 
fact that the frame tension is also a shear modulus. It 
represents the response of the system to "pure shear" de- 



formations (fig. [2fc)), which is the same as the response 
to "simple shear" (fig. Hfb)). The modulus of pure shear 
can be derived by considering the free energy cost as- 
sociated with changing the dimensions of a system from 
L x xL v x L z to [L x (l + e)] x [L y (l + e)] x [L z /(1 + e) 2 ]. 
The expansion of the free energy density in powers of e 
reads 



F 



L x L y L z 



= f a + {2P z -P x -P y )e + 0(e 2 ), (20) 



where Pj is the pressure along the i-th Cartesian axis. 
Due to the deformation, the projected area A p = L x L y 
changes by AA p ~ 2A p e. Therefore, Eq. (|2"0|) can also be 
written as 



F = F a + L Z AA P 



Pz 



P,+Py 



+ o[(AA p ) 2 ), (21) 



and the frame tension is given by 



J 



t = hm — — — = L z 

AA p ->0 AAp 



P z - 



P X +Py 



— Lz [Pn 



(22) 



r 



where P n and Pt are, respectively, the normal and trans- 
verse components of the pressure tensor (the negative 
Cauchy stress tensor) relative to the plane of the mem- 
brane. Eq. (|2"2")l is known as the mechanical definition of 
the surface tension 



25] 



D. The approximated Hamiltonian 

At this point only one question is left: why does one 
get r — do rather than r — r when dealing with the 
quadratic approximation of the effective surface Hamil- 
tonian, Eq. (|3]) ? The answer is simple - the approximated 
Hamiltonian is not rotationally invariant. This striking 



fact (which has been discussed by Grinstein and Pel- 
covitz in the context of lamellar liquid crystalline phases 
i31]) can be demonstrated by considering a certain con- 
figuration parametrized by the height function h(x,y) 
which satisfies periodic boundary conditions, and evalu- 
ating the elastic energy cost corresponding to simple and 
pure shear deformations. As discussed above, the result 
that r — fi is Hamiltonian-independent, but the following 
conclusion that fj, = t depends on rotational invariance, 
i.e., on the fact that the two shear deformations result 
in the same elastic response. This, unfortunately, is not 
the case with Eq.Q. For the simple shear deformation 
h p (x,y) = h{x,y) + ex, and upon substituting h ll (x,y) 
in Eq.Q one gets 



J 



(K {x,y))=a A p 



<l-"l.'i y(VM 2 + ^(V\ 



= nl I (h(x,y))+a A p -. 



(23) 



r 



From this result one finds that 

1 cP(n^ (hpfay))) 



A, 



de 2 



(24) 



and, thus, r = /J, = Cq. This conclusion that r = c 
is in agreement with what the equipartition theorem 



Eq.([5]) predicts for the approximated Hamiltonian. The 
response to pure shear is determined by considering the 
transformation h T (x' , y') = h(x(l + e), y(l + e))/(l + e) 2 , 
with < x' < L x (l + e), 
Ap(e) = A p (e = 
= Vft/(l + £) 3 



deformation: 
e) 2 dxdy, Vh 



<y> <L 
= 0).(l + e) 
and V 2 h T 



y (l + e). For this 
2 dx'dy' = (1 + 
V 2 /i/(l + e) 4 - 
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When these relations are used in Eq.Q, we gets: 



«»' {h r (x,y)) = (1 + t fnA p (, = 0) + / dx'dy' (Vft r ) 2 + y (V 2 ft r ) J 

<(„(*,»)) + W,<«-0)-/^[^<V*r* + **C*y|+ e *'> 
From this result, one derives the mechanical tension, which is given by 

1 d{%f{h T {x,y))) , / / , , , 



2A n 



de 



1 

-4,i 



dxdy 



<T (vhy 



k (v 2 h) 2 



r 



(25) 



(26) 



The second term on the right hand side of Eq. ([26|) is the 
entropic part of the mechanical tension, evaluated within 
the framework of the approximated quadratic Hamilto- 
nian. This entropic contribution to r is not reflected in 
the fluctuation spectrum of the quadratic Hamiltonian, 
which has the non-physical feature of not being rotation- 
ally invariant. 



IV. COMPUTER SIMULATIONS 



of the one-dimensional (ID) analogs of Eqs.© and (|3|). 
Within these two models, the membrane is represented 
by a string of N = 1024 points, the positions of which 
in 2D space are given by {fi = (xi,hi)}. In both cases, 
the simulations are performed at a constant projected 
length L p with periodic boundary conditions. Denoting 

by bi = fi + \ — fi, the distance vector ("bond") between 
adjacent points, the ID analog of the rotationally invari- 
ant Helfrich effective Hamiltonian is given by [32j 



To summarize our discussion: For any Hamiltonian 
whose corresponding free energy is given by the rota- 
tionally invariant free energy Eq.©, the correct result 
is r = t. For the not rotationally invariant quadratic 
Hamiltonian Eq.©, r = (Jq ^ t. In order to test 
these predictions, we performed Monte Carlo simulations 



J 



\bi\\bi 



(27) 



The quadratic approximation of this Hamiltonian is given 
by 



H2 — <JoLp 



op 
21 



p 4. 



i+i 



Kg 
211 



i+i 



2h l f 



(28) 



r 



where l p = Lp/N. Notice that in the ID models, the 
tension has units of a force (energy per unit length). A 
major difference between the two ID models is related to 
the positions of the points. In the rotationally invariant 
case Eq. ([2T|) . the points are allowed to be anywhere in 
the available 2D space. At low tensions, this creates con- 
figurations in which the chain forms overhangs (see e.g., 
fig. EJa)). The existence of such configurations does not 
invalidate any of the above discussion which only requires 
that the mean height profile h{x, y) is a well defined func- 
tion. Within the quadratic approximation Eq. (|28[) . the 
points are allowed to move only in the direction normal 
to the projected length [in the spirit of Eq.© in which 
the height is measured from the x — y plane]. Thus, the 
position of the z-th point is given by fi — (il p , hi) and, 



obviously, such moves do not generate any overhangs (see 
a typical configuration in fig.[2Ib)). 

In the simulations, we vary Co and measure both the 
mechanical tension, r, and the (^-coefficient, r. For each 
value of (To, the simulations extended over 2 — 4 x 10 8 
MC time units, where each time unit consists of N single 
particle move attempts and one collective "mode excita- 
tion Monte Carlo" (MEMC) move that accelerates the 
very slow relaxation dynamics of the 10 largest Fourier 
modes 33] . The introduction of MEMC moves is essen- 
tial for the equilibration of the system. The mechanical 
tension is calculated using the ID equivalent of Eq. (f22j) , 
t = f n — ft, where /„ and /* are the normal and trans- 
verse forces acting of the chain. Expressing these forces 
as thermal averages, one arrives to the following virial 
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formulae. For the rotationally invariant Helfrich Hamil- tonian "Hi 
I 



E 



(70- 



(x, + i - Xj) 2 - (hj+x - hj) 2 
l P \h\ 



Kq- 



2 / (hi - hi^i) (h i+ i - h{) (xi - Xi-if - (xi - Xi-i) (x i+ i - Xi) (hi - h^i) 2 



l P \h-i\\bi\ \ 
(xj - Xj-i) 2 + (hj - h t -iY 
N 2 

For the approximated quadratic Hamiltonian H2 

r = o*-{Y. 



\h-i\ 2 



^ (h i+ i - hi) 2 + ^ (h i+1 + hi-i - 2hi 



(29) 



(30) 



The ^-coefficient r is obtained by taking the Fourier 
transform of the the function hi, and fitting the measure- 
ments to Eq.©. Our results for r are based on the anal- 
ysis of the fluctuations of the 10 largest Fourier modes. 

Our simulation results for the quadratic Hamiltonian 
%2 are summarized in fig. [4] For both kq — and 
kq = \ksT and for all values of <7o, we find that, indeed, 
all our measurements of the g 2 -coefficient agree with the 
predicted relationship r — ao (which is denoted by the 
dashed line). As expected from Eq. (|30|) . the mechanical 
tension r is smaller than r and for small values of <7q even 



gets negative values. In comparison, the simulation re- 
sults for the rotationally invariant Helfrich Hamiltonian 
Hi are shown in fig[S] In agreement with our expecta- 
tions for this case, the various tensions satisfy the rela- 
tionship that r = t 7^ <7q. To further demonstrate the 
validity of our discussion in section Hill we computed the 
simple shear modulus /i, which for Hamiltonian T-Li, is 
also expected to be equal to r and r. The shear modu- 
lus fi can be computed using a rather cumbersome virial 
expression. For kq = 0, the virial expression for /j, takes 
the more simple form: 




Our results for /i agree perfectly with the results for r 
and r. 



V. CONCLUSIONS 

Motivated by the still ongoing debate about the various 
meanings of the term "surface tension" in membranes, we 
presented here a detailed discussion which highlights the 
importance of distinguishing between Helfrich free energy 
and Helfrich Hamiltonian, and between the Hamiltonian 
and its approximated quadratic form. Our key findings 
are the followings: 

1. We have demonstrated that, contrary to common 
perception, the observation that the spectrum of 



thermal fluctuations follows Eq.([6]), is not an evi- 
dence for the validity and accuracy of the quadratic 
Helfrich Hamiltonian Eq.©. Instead, we have de- 
rived Eq. © based on the assumption that the ther- 
modynamic properties of the membrane are cor- 
rectly depicted by Heflrich free energy Eq.©. 

2. In most thermodynamic theories the state of the 
membrane is specified by the extensive variables 
A p and A (or N) 34]. In our approach, Eq.© 
is understood as the elastic free energy functional 
that depends on the mean profile around which the 
membrane fluctuates. Thus, the state of the mem- 
brane should be described by not only two vari- 
ables, but through a function h that serves as the 
strain field for the deformed membrane. 




FIG. 3: Typical configurations of the chain in (a) the rota- 
tionally invariant Hamiltonian Eq. (|27[) . and (b) its quadratic 
approximation Eq. (|28p . Both configurations have been ob- 
tained for <7o = 0.75kBT/l p and fto = 0. In (a) the points oi 
the chain are allowed be anywhere in the available 2D space 
which creates overhangs such as those shown inside the bole 
circles. In (b) the points move only vertically to the line thai 
connects the end points of the chain. 



3. The ^-coefficient is equal to the renormalized sur- 
face tension a, i.e., the coefficient of proportionality 
between the elastic free energy and the total area of 
the mean profile in Eq.®. The renormalized ten- 
sion is a thermodynamic quantity that also depends 
on the entropy of the membrane. 

4. The result that r = a follows from linear response 
theory and, therefore, is independent of the partic- 
ular form the membrane Hamiltonian. In contrast, 
the conclusion that r = r is based on the assump- 
tion that the system is rotationally invariant and 
respond equally to all shear deformations irrespec- 
tive of the relative orientations of the deformed and 
undeformed membranes. 

5. The last point (rotational invariance) is not satis- 
fied when the membrane is described by the ap- 
proximated quadratic form Eq.Q, which calls for 
a reconsideration of some of the results derived 
through this model. Specifically, in refs. [2ll.l23.l28j 
the quadratic form has been used for a deriva- 
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FIG. 4: Simulation results for Hamiltonian H2 [Eq.p8|]: 
The ^-coefficient r (circles) and the mechanical tension r 
(squares) as a function of the intrinsic tension <to- Results 
are shown for both kq = (solid symbols) and no = lksT 
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FIG. 5: Simulation results for Hamiltonian Hi [Eq.p?])]: The 
£/ 2 -coefficient r (solid symbols) and the mechanical tension 
r (open symbols) as a function of the intrinsic tension ctq. 
Results are shown for both kq = (circles) and aco = lksT 
(squares). The pluses denote our results for the shear modulus 
\i for kq — 0. 



tion of negative mechanical tension for positive in- 
trinsic tension [see also Eq. (j2l)|) in the present pa- 
per]. Our computational results demonstrate that 
this result is achieved only with the faulty (non- 
physical) quadratic Hamiltonian H2 (see fig-H]); but 
for the corresponding rotationally invariant Hamil- 
tonian Hi, the mechanical tension is always pos- 
itive for (To > (see fig. Our computational 
results actually suggest that uo and r vanish si- 
multaneously. 

6. To express the last point in a different manner - in 
order to correctly describe the elastic properties of 
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membranes, one needs to include anharmonic terms 
in the Hamiltonian, which leads to mode coupling. 
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